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1 .  Introduction 


Mathematical  models  of  chemical  reactions  have  been  described  by  Gardiner 
McNeil,  Walls  and  Matheson  [12],  Haken  [15],  Nicolis  and  Prigogine  [ 3 1 ] , 
van  Kampen  [20],  and  Arnold  [  1  ].  For  the  deterministic  theory  of  reac¬ 
tion-diffusion  equations  we  refer  to  Smoller  [33]  and  the  references 
therein. 

In  [  3 ]  Arnold  and  Theodosopulu  have  constructed  a  space-time  jump  Markov 

process  N  by  dividing  a  finite  interval  I  (one-dimensional  reactor) 

into  N  cells,  counting  the  number  of  particles  in  each  cell  and  dividing 

this  number  by  a  proportionality  factor  v  (the  cell  size  of  an  unsealed 

model) .  This  density  changes  in  each  cell  due  to  reaction  and  diffusion 

(which  couples  neighbouring  cells) .  The  rates  by  which  this  density 

changes  are  derived  from  an  underlying  partial  differential  equation 

N2 

(PDE) .  Under  a  high  density  assumption  (—  -»  0)  Arnold  and  Theodoso¬ 
pulu  (loc.cit.)  derived  the  law  of  large  numbers  (LNN)  in  L2(I),  i.e. 

X  „  -*  X  in  La(I),  where  X  is  the  solution  of  the  PDE.  In  Kotelenez  [22] 

V,  N 

[25]  the  corresponding  central  limit  theorem  (CLT)  was  proved  under  the 
assumption  that  the  reaction  is  linear.  In  this  linear  case  the  density 
could  be  taken  low,  because  the  LLN  was  proved  in  distribution  spaces 
(cf.  also  Kotelenez  [26]).  On  the  other  hand,  nonlinear  operations  like 
multiplication  are  not  defined  on  distributions  (cf.  Schwartz  [32]). 
Therefore  it  seems  to  be  convenient  -  if  not  necessary  -  to  prove  for 
nonlinear  chemical  reactions  with  diffusion  the  LLN  in  a  function  space 
by  making  a  high  density  assumption  (as  in  Arnold  als  Theodosopulu  [  3 ]) 
and  then  derive  the  CLT  in  a  distribution  space.  This , however,  causes 
certain  numerical  difficulties  (cf.  our  Remark  3.1)  which  do  not  show 
up  if  we  assume  that  the  deterministic  limit  X  is  spatially  homogeneous 
(cf.  (2.1)  and  (2.5)).  This  assumption  allows  us  to  derive  the  LLN 
(Theorem  3.1)  in  a  function  norm  and  the  CLT  (Theorem  3.3)  in  a  distri¬ 
bution  norm.  The  limit  Y  and  the  CLT  is  a  generalized  Omstein-Uhlen- 
beck  process  (if  Yo  is  Gaussian)  and  cam  be  represented  as  the  mild 
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solution  of  a  certain  stochastic  partial  differential  equation  (SPDE) . 

We  describe  the  optimal  (smoothest)  state  spaces  for  Y.  Our  main  tool 
is  the  calculus  of  stochastic  evolution  equations  as  developed  in  Kote- 
lenez  ([2l],  [22],  [24]  -  [27])  both  for  a  fixed  Hilbert  state 
space  and  a  nuclear  Gel'fand  triple  (cf.  (2.3))  as  state  space. 

Apart  from  various  Gaussian  approximations  to  systems  of  (branching) 
Brownian  motions  (s.  our  references  in  Remark  2.2  -  and  also  Kotelenez 
[27])  we  would  like  to  mention  the  diffusion  approximations  to  spatially 
distributed  neurons  given  in  Walsh  [36]  and  Kallianpur  and  Wolpert  [19], 
where  the  limit  is  also  a  generalized  Ornstein-Uhlenbeck  process,  which 
can  be  interpreted  as  the  solution  of  a  linear  SPDE  (as  in  our  case) . 

let  us  briefly  describe  the  contents.  In  Section  2  we  introduce  both  the 
deterministic  cind  the  stochastic  models  on  an  n-dimensional  unit  cube. 

In  the  first  part  on  the  deterministic  model  we  introduce  the  nuclear 
Gel'fand  triple  (2.3)  and  prove  that  the  linear  operators  from  our  models 
can  be  "nicely"  defined  on  the  Hilbert  distribution  spaces  in  (2.3).  In 
the  second  part  on  the  stochastic  model  we  derive  some  bounds  on  N 

and  its  martingale  part.  In  Section  3  we  prove  the  LLN  in  sup-norm 
(Theorem  3.1)  with  a  certain  speed  of  convergence.  Then  we  describe  the 
limiting  Gaussian  martingale  part  for  the  normalized  martingale  parts 
of  Xv  N,  prove  in  several  steps  the  CLT  and  describe  the  limit  (Theorem 
3.3)  . 

2.  The  Models 


Following  Arnold  and  Theodosopulu  [  3  ]  and  Arnold  [  1  ]  we  first  intro¬ 
duce  the  (local)  deterministic  model,  then  construct  the  corresponding 
(local)  stochastic  model,  and  finally  compare  the  two  models. 

2.1  The  (local)  deterministic  model 


Set  S :  =  {  q  =  (q^,...,^)  €  IR^  :  0  £  q  S  1 ,  i  =  1 , 


•  •  •  » 


n)  .  Let 


R(x)  =  b(x)  -  d(x)  =  c^x  be  a  polynomial  in  x  €  3R,  where  2  0, 

c  <0  and  b(x)  and  d(x)  are  polynomials  of  degree  <  m  with  nonnegative 
m  - 

coefficients.  A  denotes  the  Lapacian  and  D  >  0  a  diffusion  coefficient 
Then  the  concentration  of  one  reactant  with  reflection  at  the  boundary 
is  given  by  the  solution  of  the  following  PDE: 


—  X(t,q)  =  oAx(t,q)  +  R(X(t,q) ) 


(2.1)  c  —  X(t,q)  *0,  q.  6  {0,1}  i  =  l,...,n 
dq .  1 


Xo(q)  >  0 


Let  Bq  :=  L2(S)  be  the  Hilbert  space  of  square  integrable  real  valued 

functions  on  S  equipped  with  the  scalar  product  <cp,iji>0:=  fcp(q) ii»(q)  dq, 

s 

cprxjd  6  E  .  In  what  follows  we  shall  denote  by  DA  the  closure  of  DA 
o 

w.r.t.  the  reflecting  boundary  conditions  of  (2.1).  DA  is  self-adjoint 
nonpositive  on  Eq  and  has  a  discrete  spectrum.  Let  l  =  (£^,...,£  )  be 

a  multiindex,  where  £  6  IN  U  {o} ,  and  set 


^2  cos  JLnf*) 


l.  2  1 
1 

l.  =  0 

X 


Then,  the  :=  ^11^  <t> ^  are  a  complete  orthonormal  system  (CONS)  of 

1  £}  2  2 

eigenvectors  of  DA  with  eigenvalues  -DU^  :  =  -D(  Zj  JL  Tt  )  .  Consequently 
the  semigroup  T(t)  generated  by  DA  on can  be  represented  by 


(2.2)  T(t)<p  =  Z  e 
l 


"D%t 


l  ' 


As  in  Kotelenez  ([26],  [27])  we  introduce  the  nuclear  Gel' f and  triple 
determined  by  DA 


(2.3)  ®  c  E  c  H  =  B'  c  H  _  C  ©•  ,  a  >  0. 


In  (2.3)  we  have  : =  Dom ( (I  -  DA)  ) ,  a  £  0  where  I  is  the  identity 

operator  and  "Dom"  denotes  "domain".  B^  is  a  real  separable  Hilbert 
space  if  equipped  with  the  scalar  product  <♦,•>  :=  <(I-DA)<3^^* ,  (I-DA)01^2* 

(for  the  definition  of  the  a-th  power  of  a  positive  self-adjoint  opera¬ 
tor  -  cf.  Yosida  [38]).  E’,  the  strong  dual  of  E  ,  is  identified  with 

o  o 

E  ,  $  =  fl  K  is  a  locally  convex  vector  space  whose  topology  is  given 
°  c£0  a 

i 

by  the  set  of  norms  {Icpl^  :=  (<cp,cp^  )^"  ,  cp  €  ®]  ,  and  is  the  strong 

dual  of  ®.  E  are  those  <p'  £  ®  which  can  be  extended  to  continuous 
-CL 

functionals  on  E^  ,  a  ^  O.  E^  is  3  real  separable  Hilbert  space 
with  scalar  product  <• , •>  ,  where  for  CP, 

-CL 

cp  £  E  <  cp,ip  >  =  <  (I-DA)-Cl^2<p,  (I-DA)-0^2^  >  .  Moreover,  setting 

o  -cl  o 

:=  1  +  DU^  ,  we  obtain  that 


is  a  CONS  for  E  ,  a  €  E.  Hence 
a 

Ba  =  {cp1  £  :  I  (CP’,$£)2  A.®  <  oc  }  , 

where  ( • ,  • )  denotes  the  dual  pairing.  Thus,  if  we  set 

Yu  -  “v  e  m~  >  f  4  4  <-> 

oo 

we  see  that  (2.1)  can  be  identified  with  a  subset  of  E  ,  where  E  is 

a 

isomorphic  to  1  ,  a  £  E.  Clearly,  the  imbeddings  in  (2.3)  are  con- 

Z 

tinuous  and  dense . 

Lemma  2 .  1 

For  any  a,y  £  E  s.t.  a  >  y  +  -j  the  imbedding 


is  Hilbert-Schmidt. 


ov 


Proof 


i  !♦;  i’  =  i  i*a-v  <-  ^ 

oo  oo 

f  ...  f  (1+x2  +  ...+  X2)^+Y  dx.  . . .  dx  <  oo  iff 
^  '  1  n  In 

oo  oo 

/ _ /  (l+x1  —  +  x^)  2a+2'f  dx^  ...  dxn  <  <*>  iff 

2a  >  n  +  2y. 

a 

a/2 

Since  (I-DA)  and  T(t)  commute ,T (t)  can  be  extended  (resp.  restricted) 
to  a  strongly  continuous  semigroup  T^(t)  on  E^,  a  €  E,  s.t.  for  all 
a  €  E 


(2.4) 


It  (t) I 


L  (H  ) 

a 


It ( t) I 


L  (H  ) 
o 


<  1. 


L(Ha)  denotes  the  usual  operator  norm  on  BQ,  and  the  inequality  in 

(2.4)  holds  because  DA  is  dissipative  (cf.  Davies  [9]).  Let  us  denote 
by  DA^  the  generator  of  T^(t)  (which  is  the  extension  (resp.  restriction) 

of  DA).  As  in  Kotelenez  [26],  Lemma  2.2,  we  obtain: 


Lemma  2 ■ 2 

For  all  a  €  E,  Dom  (DA  )  =  E  _  and  T  (t)  is  analytic. 

a  a+2  a  1 

For  the  rest  of  the  paper  we  shall  assume  that  the  system  (2.1)  starts 

in  a  spatially  homogeneous  state  XQ  =  pD  >  0.  This  implies  that  the  solu¬ 
tion  X(t,r)  of  (2.1)  is  spatially  homogeneous,  i.e.  X(t,r)  e  p(t)  satis¬ 
fies  the  ordinary  differential  equation 

(2.5)  |rP(t)  =  R(p(t) ) ,  p(0)  «  P  >  0 

dt  o 

(cf.  Arnold  [  1  ]).  (2.5)  has  a  unique  positive  bounded  solution  which 
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is  strictly  positive  for  all  t  >  0  (cf.  Coddington  and  Levinson  t  5 ] ) - 

2.2  The  (local)  stochastic  model 

We  cover  S  with  grid  of  N  n-dimensional  cubes  (cells)  of  size  hn  which 
are  parallel  to  the  axes.  The  cell  corresponding  to  the  grid  point  r ^  is 
defined  by 


[ r J )  :  =  {r  €  S:  r"  <  r.  <  rJ  +  h,  i  =  1,  . 

>  ill 


•  * i  nj  t  ^  If  « . «  /  N . 


Let  v  be  a  parameter  (which  is  explained  in  Remark  2.1)  and  denote  by 

TE  the  (countable)  state  space  of  elements  k  =  (k  -J.r  where 

N  rJ  ( L  rJ  )  cs ) 

k  ,  €  —  U  {O}  .  Set 
rJ  v 


E  . ,  {cp  €  E  :  cp  constant  on  each  [r^)  }  . 
O ,  N  O 


r: 


E  cz  ffl 
N  o,N 


Tt  :=  E  -»  E  defined  by 
N  o  o,N 


ti  <p(r)  :  =  cp  (r )  :=  h  n  /  cp  :=  h  n  J.  cp(q)dq  if  r  £  [r^) 

N  N  [rj>  Er^) 

is  a  projection  from  E  onto  E  . 

o  o,N 

Now  we  define  a  Markov  chain  on  3E..  through  the  Q-matrix  of  its  transi- 

N 

tion  intensities: 


r  3 5 

m 

=  k 

+  er: 

=  :m 

+o 

ri) 

m 

=  k 

"  V 

=  :  m 

-o 

,0 

kr3 

m 

=  k 

+  e  -e 

r3±hi  rj 

= :  m  , 
±i 

(2.6)  B(k,m)  :=  ' 


-.1  &(k,m±i) 

i=o  11 


for  O  <  ±  h.  ^  1 

l 

m  =  k 
otherwise . 


‘J 

.i 


Here  e^j  =  d[rj))  *  where  l[rjj  (q)  =  1  if  q  €  [r3),  =  O  otherwise  , 

and  h^  =  (0, . . ,0, h,0, . . ,0)  where  all  but  the  i-th  coordinate  are  zero. 

Hence,  we  obtain  the  distributions  P(t,k)  determined  by  Q  =  (3(k,m)  _ 

k  ,m€  IE 

as  the  unique  solution  of  Kolmogorov's  backward  equation  (which 
is  called  in  the  application-oriented  literature  the  "multivariate 
Master  equation)  (cf.  Arnold  [  1  ]).  The  corresponding  (canonical)  cadlag 
Markov  process  will  be  denoted  by 

(generated  by  Q) 

(2.7) 

X  .,(0)  =  X  £  E  ,  a  given  initial  distribution, 

v ,  N  v,N,o  N 

In  what  follows  we  shall  assume  that  the  stochastic  basis  for  X  „ 

v,N 

(fi,  F,  F  ,P)  is  complete  with  right  continuous  filtration. 

Remark  2 .  1 

We  can  view  Xv  N  as  the  rescaled  density  Markov  process  of  Arnold  [  1  ] 
and  Arnold  and  Theododopulu  [  3 ]  on  a  cube  of  volume  vN  =  V  with  cells 
of  size  v,  where  the  number  of  particles  is  proportionel  to  v. 

Remark  2 , 2 

If  b(r)  =  br  +  c  ,  d(r)  =  dr,  for  some  constants  b,d  >  0  then  X  „  is  a 
branching  diffusion  with  immigration  (c^)  on  the  grid.  This  case  was 
investigated  in  Kotelenez  (25],  [26],  and  the  limit  theorems  therein 
corresponded  to  limit  theorems  for  branching  Brownian  motions  obtained 
by  Holley  and  Stroock  [16],  Gorostiza  [ 1 4 ]  and,  in  the  absence  of  branch¬ 
ing  (b  =  d  =  c  =  0)  to  Martin-Lof  [30]  and  Ito  [18]  (cf.  also  Walsh 
o 

[3"]).  For  a  diffusion  approximation  to  branching  Brownian  motions  -  cf. 
Dawson  [11]. 


In  what  follows  we  shall  not  explicitly  write  the  parameter  v,  i.e.,  we 

shall  write  X,,  instead  of  X  etc. 

N  v,  N 


Extend  cp  ,  £  E  by  reflection  to 
N  o ,  N 


:=  {r  £  E  :  -h  <  r.  <  1+h,  i  =  1,.. 
h  i 


and  set 


V^i  «>N  :=  h_1[(PN(r±h.)  -  cpN(r)] 


A  CP  -.£,***  V%  . 

N  N  1=1  -i  i 


Remark  2 . 3 


In  view  of  our  boundary  condition  we  easily  see  that  DA^  is  selfadjoint 
and  dissipative  both  as  an  operator  on  Eq  n  and  Eq  (where  on  A^  is 

defined  by  A„oTT  )  .  If  we  set 
N  N 


l.  ,N  ° 

l 


we  see  that 


{tPIfN  :=  iSl  V  ,N  '  *i  <  h  1 


is  a  CONS  of  eigenvectors  of  DA  for  E  with  eigenvalues 

N  o,N 


DU.  :=  2nDN^.n,  {l  -  cos  i .  hTt}  . 
£,N  i=l  i 


The  waiting  time  parameter  for  XN  is  given  by 

a  (k)  =  v  I  l R I (k  j)  +  h  2.I.  2Dk  j 
N  r:£s  r  i-l 

with  IrI (x)  =  b(x)  +  d(x) ,  x  €  1R.  Hence,  if 


0..(k,m)  :=  (0..(k))  B(k,m) 


denotes  the  jump  distribution  function  (&(k,m)  from  (2.6)),  then  the  in¬ 


finitesimal  generator  for  X„  is  given  by 

N 


(2.8)  (Af) (k)  =  Of k)  /  [f (m)  -  f(k)J  0M(k,dm) 

N  *  N 


N 


where  f:  E  -*  3R  is  bounded  and  measurable  (Gihman  and  Skorohod  [13]). 
N 


Let 


<J>  Ill  :*  sup  l<(>  (r)  I 
N  r€S  N 


be  the  sup-norm  on  E  ..  If  there  is  a  finite  constant  K(v,N)  s.t. 
*  0,N 


(2.9) 


III  X  „(OJ  III  <  K (v, N)  a.s. 
v,N 


then  by  a  lemma  of  Kurtz  [28]  (cf.  Arnold  and  Theodosopulu  [  3  ]  and 
Kotelenez  [25]) 


'yt}  ==  yt>  -  xN(0)  -  /  aN(xN(s))J  (z-  xN(s))0N(xN(s),dz) 


ds 


E, 


N 


(2.10)  *!  t 

=  X  ( t )  -  X  (0)  -  J  [DA  XM(S)  +  R(X  (s))]ds 
N  N  N  N  N 

O 


is  an  H  -valued  square  integrable  cadlag  martingale. 

o,N 


We  shall  assume  (2.9)  throughout  the  paper. 


Hence,  X  satisfies  formally  the  stochastic  evolution  eqaution 
N 


f dX  (t)  =  [DA  X,(t)  +  R(X  (t) ) ]dt  +  d2  (t) 
N  N  N  N  N 


(2.11) 


x.»  ( ° )  =  xM 

N  N,o 


and  the  difference  X  (t)  -  X ( t) ,  where  X  is  the  solution  of  (2.1)/ (2.5) 

N 


satisfies 


-  10 


fXN(t) 


(2.12) 


-  X(t)  =  XN(0)  -  X(0)  +  /  (DA^  R'  (X(s))  (XN(s)  -  X(s)  )ds 


+  /‘'(XN(s)  -  X(s)  ) 2  R(Xn(s)  ,X(s))ds 


+  2N(t) 


R' (x)  is  the  derivative  of  R(x) ,  R(y,x)  is  a  polynomial  in  y  and  x  of 

2 

degree  £  m-2,  and  R'(X(s))  and  (XN(s)  -  X(s))  are  interpreted  as  multi¬ 
plication  operators. 


Note  that  both  DAj^  +  R'(X(s))  and  DA  +  R'(X(s))  are  quasi- generators  of 
evolution  operators  UN(t,s)  and  U(t,s)  on  ^  and  E^,  respectively. 
(For  the  definition  of  evolution  operators  V(t,s),  i.e.,  strongly  con¬ 
tinuous  two-parameter  semigroups  -  cf.  Curtain  and  Pritchard  [  7  ]  and 
Tanabe  [34],  where  V(t,s)  is  called  fundamental  solution  -  any  strongly 
continuous  one-parameter  semigroup  is,  of  course,  also  an  evolution 
operator.)  Consequently,  by  variation  of  constants,  (2.12)  yields 


vt}  - 


X(t)  =  uN(t,0) (XN(0)  - 


rt 

X(0))  +  J  U  (t,s)dZ  (s) 
o  N  N 


(2.13) 


+ 


/CUN(t,s) (XN(s)  -  X(s))2  R(X^(s) ,X(s) )ds. 
o 


In  order  to  give  a  meaning  to  the  stochastic  convolution  integral  in  (2.13) 
we  recall  from  Kotelenez  [2l],  [24]: 


Definition  2.1 

Let  H  be  a  separable  Hilbert-space  with  Hilbert  space  norm  I*  I  and 

H 

V(t,s)  an  evolution  operator  onH,  0  £  s  ^  t  <  <*>.  V(t,s)  is  of  contrac¬ 
tion-type  or,  equivalently,  V(t,s)  6  G(l,&*)  if  for  all  t  >  0  there  is 
a  finite  constant  >  0  s.t. 


(t-s) 


(2.14) 


IV(t,s)  I  ,  ,  <  e 


11 


•  '  -  -  •'  •  ■'*  -  *-T  -.-V  -  n.’.  V."  K"  "‘.I KT.  - 


for  all  0  <  s  £  t  £  t. 

Remark  2.4 


Let  M  be  an  E-valued  locally  square  integrable  cadlag  martingale  and 
V(t,s)  £  G(l,3£)  onE.  Then,  from  Kotelenez  [21],  we  have 

(i)  If  M  is  cadlag,  then  J  V(*,s)dM(s)  has  a  cadlag  version?  if 

o 

M  is  continuous,  then  /  V(*,s)dM(s)  has  a  continuous  version. 

A  partial  result  from  Kotelenez  [24]  is  the  following: 


(ii)  If  V(t,s)  has  a  quasi-generator  A<t)  and  Dom(A(t))  is  independ¬ 
ent  of  t  then  for  all  t  >  0  there  is  a  finite  constant  c  =  c(t,3) 

depending  only  on  the  scalar  product  <  *,*  >  ,t  and  3  s.t.  for 

H 

all  t  ^  t 


(2.15) 


E  sup  I J  v(s,u)dM(u) 1^ 
o<sSt  o  H 


< 


c  e 


43-t 


ElM(t) I 


2 

H  • 


For  more  general  properties  and  inequalities  for  stochastic  convolution 
integrals  cf.  Kotelenez  [24],  [27]. 

Since  X (t)  is  constant  in  the  space  variable  (cf.  (2.5))  we  obtain 

f  UN(t,s)  =  (t-s )  exp ( J  R'(X(u))du) 

(2.16)  <•' 

U  ( t , s )  =  T  (t-s)  exp(/  R ' ( X ( u ) ) du ) ,  a  €  2R, 

—  s 

where  the  last  equation  means  that  U(t,s)  is  extendible  (resp.  restrict- 

able)  to  the  E  .  Let  us  denote  by  (  E  .<•,•>)  E  „  equipped 
a  J  a,N  a  o,n 


with  the  Hilbert  norm  <•••>,  a  £  1R.  Set 


‘.Vi 


12 


0  ;  =  sup  R' (X(t) ) 
o£t<® 

and  note  that  0  <  <*>  by  our  assumption  on  R(x)  (cf.  (2.5)). 


Lemma  2 . 3 


For  all  a  €  3R 


U  (t,s)  6  G(1,0)  on  Ea 


V  At, s)  €  G(1,0)  on  E  M 

N  OL»N 


Proof 


(i)  The  statement  for  U^t/S)  follows  from  (2.4)  and  (2.16). 


Let  x  €  E  „  .  Then 
o,N 


luN(t's)xla  =  |  <  Vt,s)x'*Jl  >2o 


=  I  <  x,uN(t,s)VN  >oX? 


<  ie2&{t-S>  <  x,TN(t-s)^jN  >2o  X® 


by  (2.16) 


since  T  (t)  is  a  contraction  and  „  is  an  eigenvector  of  T  (t) • 

N  x,  ,N  N 


The  previous  considerations  show  that  the  limit  behaviour  of  X^(t)  "  ^(t) 
essentially  depends  on  the  limit  behaviour  of  Z^(t),  U^(t,s)  and  the 


last  term  in  (2.13). 


We  shall  first  give  an  estimate  on  the  variance  of  2  (t) .  To  this  end  we 

N 

define  an  operator  on  by 

(2.17)  F  (tp)  :  =  D  .1,  V*1.  +  V*  Cp J*.  +  |R|(cp) 

N  1=1  -li  l-l 

where  cp  €  and  i R I  (cp)  act  as  multiplication  operators.  As  in  the 
linear  case  of  Kotelenez  [25]  we  obtain 

Lemma  2.4 

For  arbitrary  cp  £  H 

o 

(2.18)  *  <  VCI '®  >o  *  h  £  E<F»<XKl5>>V»N>oas- 

We  need  estimates  on  XN(t,r)  which  satisfies  by  variation  of  constants 

(2.19)  XN(t)  =  TN(t)XN(0)  +  ^  TN(t-s)dZN(s)  +  J*  TN(t-s)R(XN(s))ds 

Set  p  :=  max  {R(x)  :  x  €  JR+]  . 

The  definition  of  R(x)  implies  p  <  =>. 

Lemma  2 . 5 
For  any  t  >  0 

(2.20)  sup  III  E  X  { s )  III  £  tp  +  III  E  X  (0)  III  , 

oSsit  N  * 

(2.21)  III  X  (t)  III  <  tp  +  III  X  (0)  III  +  /n  I  T  ( t-s ) dZ  ( s )  I  . 

N  N  q  N  N  o 

Proof 

(i)  From  Kotelenez  [22],  lemma  A. 7  and  Davies  [9  ],  Th.  7.16  we 
obtain  that  T  (t)  is  positivity-preserving  on  E  ,  i.e., 

O  t  N 

leaves  the  cone  of  nonnegative  functions  invariant,  which  implies 


(ii)  We  easily  check  that  for  any  <p„  £  E  „ 

N  O  ,N 

(2.22)  III  III  <  &  1(0,1  , 

N  NO 

whence  we  obtain  (2.21)  from  (2.19). 

□ 


3.  Limit  Theorems 


Theorem  3 . 1  (LLN) 

Assume  (2.9)  in  addition  to 

2 

2y+  1-t-jj- 

(I)  v  =  NP,  where  P  >  — - —  and  y  6  [  j  ,  y  )  arbitrary  and  fixed 

(II)  E  III  (vN)Y(X  (0)  -  X (O) )  III  -»  0,  as  N  -*  •  . 

N 

Then  for  all  t  >  O,  6  >  O 


P{  sup  III  (vN)Y(X„(t)  -  X(t)  III  >  6}  -»  0;  as  N  -♦  oo  . 
o<t<t  N 


Proof 

(i)  (2.18),  (2.20)  and  our  assumptions  imply  the  existence  of  a  finite 

constant  K  s.t.  for  any  t  >  0 


(3.1) 


iz  (t)  r 

N  o 


K(t+1) 


(c  +  N2/n) 
o 


where  by  (2.15)  and  Lemma  2.3  there  is  for  any  t  £  0  a  finite 

~  r  1  1 
constant  K(p,t)  s.t.  for  y  £  l  —  ,  —  )  and 


nM(t)  :=  (VN)YN1/2  max  {|/fcU  (t,s)d2  (s)l  , Ij*  T  (t-s)dZ  (s) I  } 


max 


E 


-  15  - 

sup  nj(t)  <  Kte.t)  -^V'N--  N(c  +  N2/n) 
o<t<t  ™  VO 

=  KCfJ.t)^2^  +  2/n  +  1  +  2Y 
-»  O,  as  N  -*  oo'  by  the  definition  of  P- 

(ii)  (2.16)  implies  that  U  (t,s)  is  positivity  preserving  since  T„(t) 

N  N 

is  positivity-preserving.  Abbreviating 

^(t)  :=  e&t  (vN)Y  III  XN(0)  -  X(O)  III  +  ^(t) 

and 

lk,(t)  :=  e0t  III  E(X  (t),  X(t)  III  •  III  X  (t)  -  X(t)  III  ,  t  S  t 

N  N  N 

the  Gronwall-Bellmann  lemma  and  (2.13)  imply 

(vN) Y  III  xN(t)  -  x(t)  III  £  ^(t)  +  CN(s)^N(s)e3q>(/t^fJ(u)du)ds 

t  £  t. 

Since  by  step  (i)  sup  4/  (t)  is  stochastically  bounded  as  N  ■*»  and 
o£t£t  N 

sup  LT(t)  tends  to  zero  in  mean  square  the  proof  is  finished. 
o<t<t 

a 

Set  1/2 

Mn  (VN)  ZN 

and  define  for  cp  €  the  continuous  analogue  to  (2.17): 

(3.2)  F«p)  :  =  -2D . £  3. CP  3.  +  |R|  (CP) 

i-l  i  i 

where  again  cp  and  |R|  (cp)  act  as  multiplication  operators.  Denote  for 
U  €  (0,1)  by 


CU([0,co)  ;  H) 


the  space  of  HAlder  continuous  E-valued  functions  with  Haider  exponent  U 
where  E  is  some  Hilbert  space. 


Lemma  3 . 1 


There  is  a  unique  {in  distribution)  ®' -valued  Gaussian  martingale  M  on 

~  ^ 

some  probability  space  (£l,F,Ft,P)  with  characteristic  functional 

~  i 

(3.3)  E  exp  (i  (M(t)  ,<p) )  =  exp(-  -z-  J  <  F(X(s))cp,<p  >  ds) , 

*  o  ° 

_  ~  row 

ip  t  4,  where  E  is  the  mathematical  expectation  vr.r.t.  P.  Moreover,  for 
any  a  >  j  +  1  and  any  y  €  (0,  j) 

M  €  CU([0,»)  ,  E  )  a.s. 

-CL 

The  proof  of  the  existence  and  uniqueness  is  given  in  Ito  [17] (cf .  also 
Ustunel  [35]),  and  the  H61der  continuity  follows  from  Kotelenez  [23]. 

Since  X(t)  is  spatially  homogeneous  and  strictly  positive  we  easily 
check  that  F(X(t))  as  a  positive  self-adjoint  operator  on  Eq  is  just 
equal  to  -  2X(t)DA  +  lRl(X(t)),  which  has  the  same  eigenfunctions  <t>^  as 
DA.  Thus,  the  square  root  of  F(X(t))  can  be  considered  as  an  element 

1  /o 

F  (X (t) )  from  L(  B  ,  E  .)  for  all  a  €  E.  If  a  >  r  then  there  is  an 
a  a  a- 1  2 

E_a  valued  Wiener  process  W(t)  on  E_^  which  is  the  cylindrical  Brownian 

motion  on  E  (cf.  It6  [17]).  We  may  without  loss  of  generality  assume 

that  W (t)  is  also  defined  on  (ft,  rrF^,P)  from  Lemma  3.1.  Repeating  now 

the  proof  of  Lemma  2.4  in  Kotelenez  [25]  we  obtain 

Lemma  3 . 2 

(3.4)  M  =  /  (X(s)dW(s)  (equal  in  distribution) 

o 

on  CU  ([0,eo);  K  )  for  all  a  >  §■  ,  all  u  £  (0,  . 


Let  us  denote  by  D([o,<*>);  E)  the  complete  metric  space  of  E-valued  cadlag 
functions,  where  E  is  a  separable  Hilbert  space  (i.e.  the  Skorohod  space 
-  cf.  Billingsley  [  4  ]  and  Kurtz  [29])  and  by  *W  weak  convergence. 


Lemma  3 . 3 

Under  the  assumptions  of  Theorem  3.1  for  all  a  >  y  +  1 

M  •  M  on  D([0,oo);  E  ), 

N  —CL 

where  M  is  the  Gaussian  martingale  given  in  Lemma  3.1. 


Proof 

(i)  The  weak  convergence  of  M  ft)  to  Mft)  for  fixed  t  follows  as  in 

N 

the  linear  case  (cf.  Kotelenez  [25],  [26]). 

(ii)  We  shall  estimate  the  "modules  of  continuity".  Set  for  some 
(large)  K  >  0 

T  :=  inf{t  >  0  :  III  X  (t)  III  >  K} 

N  N 


and  MN(t)  :=  ^(tA^) 

where  "a"  denotes  "min".  Then,  abbreviating  F  :=  a(XM(s),sSt) 

a  N ,  t  W 

we  obtain  for  t  £  t,  s  >  0 


(*) 


: E  ii  fk>i) 

-  s  E  "V1*5’  FBjt) 


*  E  UhN  <  its)  1  FN,t’  ' 


Take  the  CONS  {<?“  }  for  E  .  By  (2.18)  (cf.  Kurtz  [28]  and  Kote¬ 
lenez  [22],  [25]  for  the  step  from  the  unconditional  to  the  con¬ 
ditional  expectation)  the  first  term  in  the  r.h.s.  of  (*)  can  be 
estimated  from  above  by 


!=(/  <  v  VV“V’*I>od“ 'Vt1 


by  Lemma  A. 2  in  Kotelenez  [22] 

i*«< 

<  Ks 

for  some  K  <  «  since  by  Lemma  2.1  E  ‘ — *  E  is  Hilbert-Schmidt 

a- 1  o 

and  III  X  ( tAT  )  III  <  K  +  1  <  oo  . 

N  N 

Setting 

Vtl!l  +  i{T„  <  L.) 

N 

we  obtain  from  Theorem  3 . 1 

lim  lim  E  y  *(s)  -  0. 
s-o  N-«o 

(iii)  (i)  and  (ii)  imply  by  Theorem  2.7  of  Kurtz  [29]  the  weak  conver¬ 
gence  of  M  to  M. 

N 


Theorem  3.2 

Under  the  assumptions  of  Theorem  3.1  for  all  a  >  y  +  1 

/  U  ( • , s ) dM  ( s )  -  j’u  (•  ,s)dM(s)  on  D([0,<»);E  ). 

N  N  -CL  H3, 

O  O 

Proof 

(i)  Let  dp  denote  the  Prohorov  metric  on  D([o,“>);  E^)  (cf.  Billingsley 

[4]).  Let  tt  be  the  projection  of  E  onto  L(<p  :  L.  <  k  for  all 
K  “CL  i 

i  *  l,..,n)  (the  linear  hull  spanned  by  those  whose  multiindices 


I  =  (£^,...,£  )  satisfy  £  <  k  for  all  i  =  The  corresponding 

projection  f  rom  E  onto  L  ( <f>  „  :  l.  <  k  for  all  i  =  l,...,n)  if 
o,n  a , N  1 

-1  N  - 1  N 

k  <  h  will  be  denoted  by  py..  If  k  >  h  then  we  set  pv  B  „  =  B  . 

1  **  k  o,N  o,N 

1  IN  N 

Set  n  *  I  -  and  p^  =  I  -  p^,  where  I  denotes  the  identity  opera¬ 
tor  on  the  corresponding  spaces. 


(ii)  Abbreviating  the  convolution  integrals  /  U  (•  ,s)dM^(s)  by  l^dM^  etc 

o 

we  obtain 

dp(/UNdMN,  /0dM) 

<  d  (fu  dM  ,  fu  pNdM  ) 

N  N  j  N^cV 

(*) 

+  V^V^V  ^UTlkdM) 

+  dp( Jun^dM,  JudM) . 


(iii)  By  Lemma  2.3  and  (2.15)  for  any  t 


E  sup_  I/1  0 
o<t£t  ° 


s  =  ,4Bt  El  A.(i) 


»k  VtM-a 


_  A 

S  c  e4Pt  I  E  <  M  (t)  ,  4>«n  >1 

a  >k  *  *N  ° 

i 


^  c  e4^  K ( t+ 1 )  I  X~a+1 
£.>k 


>  0 


as  in  the  proof  of  Lemma  3.3  for  some  finite  constant  K. Since 

I  k?*1  <  oo  by  Lemma  2.1  the  r.h.s.  of  the  last  inequality  can  be 
£  1 


made  arbitrarily  small  by  choosing  k  large.  Hence  for  given  E  >  0 
there  is  a  k. (e)  s.t.  for  all  k  £  k. (E)  and  all  N 


V^W  J’vS’V  s  ! 

(cf.  Kotelenez  [24],  [25]). 

The  third  term  in  (*)  can  be  estimated  in  the  same  way.  The  second  term 
in  (*)  tends  to  zero  for  fixed  k.  Indeed,  by  partial  integration 

/tUN(t,s)P^dMN(s)  =  p^MN(t)  +  J  UN(t,s)[DA^  +  R' (X(s) )  ]p^MN(s)ds  and 

JtU(t,s)TljtdM(s)  =  Tt^M(t)  +  /t  U(t,s)  (DA  +  R*  (X(s)  )TT^M(s)ds.  Hence,  the 
o  o 

Trotter-Kato  theorem  (Davies  [9],  Theorem  3.17,  Kotelenez  [24],  Remark 

N 

4.1,  and  Kotelenez  [22],  Lemmas  A.l,  A. 3)  and  the  definition  of  pk  and 
Tt,  imply  the  conditions  of  Theorem  5.5  in  Billingsley  [4],  Ch.I. 


(iv)  Since  weak  convergence  on  D([o,<»);  B  )  and  convergence  w.r.t.  the 


-a 


Prohorov  metric  dp  are  equivalent  (cf.  Kurtz  [29]  and  Billingsley  [  4  ], 
Appendix  III,  Th.  5)  the  proof  is  finished. 


Fix 

a  >  |  +  1  . 

~  ~  ~  ~ 

Let  Y  be  an  E  -valued  square  integrable  random  variable  on  (Q, r,r  ,P) 
o  -a+ 1  t 

independent  of  W(t)  for  all  tt  0  and  Y^  a  square  integrable  ®a+1 -valued 

random  variable  on  (£2,  F,F^,p)  such  Yq  -  Y^.  Further,  let  C  denote  an 

arbitrary  / U  F^2^  dW  -  continuity  set  of  D([o,«>);  B^)  (cf.  Billingsley 

[  4  ])  and  E  an  arbitrary  element  from  a(Y^) .  We  make  the  following 
asymptotic  independence  assumption: 

'  Y  :=  (vN)1//2(X„  -  X  )  -*  Y  in  probability  on  E 

N ,  o  N ,  O  o  o  -vl+ 1 


(cf.  Billingsley  [4],  Ch.I,  Th.  2.1  -  The  second  condition  in  (3.5)  is, 

e.g.,  satisfied  if  Mv_  is  independent  of  O(Y^)  for  all  (large)  N. 

N  u 

Let  6  denote  the  Frechet  derivative,  B([o,t]  x  E  )  the  real  valued  mea- 

-  2 

surable  functions  g  with  domain  [0,t]  *  E  s.t.  -rf-  ,  6g,6  g,  and 

“CL  at 

Dh_^6g  exist,  are  continuous  in  x  and  t,  and  uniformly  bounded  in  norm 
1/2 

on  [0,T]  x  E  .  Q  is  the  square  root  of  the  covariance  operator  of 
—CL 

1/2*  1/2 
W(t)  on  E  ,  and  F  ,  (X(t))  is  the  dual  operator  of  F  (X(t)) 

-0L+ 1  -CL+ 1  “CL+ 1 

(after  identifying  the  duals  of B  and  E _  ,  with  E_^  and  E  ,  re- 

-a  -oi+ 1  —a.  -ci+ 1 

spectively.  Finally,  "Tr"  denotes  "trace". 

Now  we  can  state  our  final  result  under  the  assumptions  of  the  LLN. 


Theorem  3.3  (CLT) 

Assume  (2.9)  and  (3.5)  for  fixed  a  >  j  +  1  in  addition  to 

2y+l+^- 

(I)  v  =*  IT  where  p  >  i-fy*  3111(3  Y  £  ^  ?'  2*  arbitrary  md  fixed 

(II)  E  I  (vN)Y(X  (0)  -  X(0))|2  -»  0,  as  N  . 

N  O 

1/2 

Then  for  YxT  :=  (vN)  (X  -X) 

N  N 

(i)  yN  *Y  on  D ( [0 ,°°)  ;  E_^  ) 

where 

(3.6)  Y(t)  =  U  (t,0)  Y  +  ft  U  ( t ,  s )  F1^2  (X(s) )  dW  (s) 

-a.  o  q  -cl  -a,+ 1 

is  the  mild  solution  of  the  stochastic  partial  differential  equation 

|  dY(t)  =  [DA  +  R’ (X(t))]Y(t)dt  +  F1^2  (X(t) ) dW (t) 

1  -a  -0L+ 1 

Y (O)  =  Y 


(3.7) 


|  Ci  [0,t] ;  E^)  a.s.  for  all  U  <  —  ,  all  t  >  0 
(3.8)  Y  €  < 

C( [o, t] ;  E  )  a.s.  for  all  t  >  O. 

-a+ 1 


and  Y ( t) ,  t  >  O,  does  not  define  a  G-additive  measure  on  E  for  y  £ 

-01 

i.e.,  the  second  relation  in  (3.8)  is  the  maximal  regularity  of  Y  on  the 
Hilbert  scale  (2.3). 


(iii)  Y  is  a  Markov  process,  and  its  weak  generator  is  given  by 


A  (t)  g  (t,cp* )  =—  g(t,<p')  -t-  <[DA_a  +  R' (X(t) )  ]6g(t,(p' )  ,cp’> 


-a 


(3.9) 


+  jTr{Q1/2  F^2*(X(t))62g(t,ip*)F^21(X(t))Q1/2}  , 


where  g  €  B < [O, t]  *  E  ) . 

-a 


Proof 


(i)  The  norm  of  the  normalized  last  term  in  (2.13)  can  be  estimated 
as  follows: 

I /fc  U  ( t ,  s ) (X  (s)  -  X(s) ) 2(vN) 1/2  R(Xm(s) ,X(s) ) I 

N  N  N  -a 

o 

<  eB(t"S)  111  (vN)  1/4(Xm(s)-X(s)  )  III2-  III  R(X„(s' ,X(s)  III  ds  -»  0, 

N  N 

O 


as  N  -*  oo  uniformly  on  compact  intervals  in  probability  by  Theorem  3.1 

and  the  stochastic  boundedness  of  III  R(X  (s)  ,X(s)  III  (cf.  the  proof  of 

N 

Theorem  3.1).  Therefore,  the  weak  convergence  of  YN  to  Y  follows  from 
Theorem  3.2  and  our  assumptions  as  in  Kotelenez  [25],  [26]. 


c-  .*■ 


NJ|D 
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(ii)  The  HSlder  continuity  follows  from  DaPrato,  Iannelli  and  Tubaro 
[  8]  and  Kotelenez  [23]. 

The  spatial  regularity  follows  from  the  estimate 

if  U^(t,u)F^f  (X(u))dW(u)l<x+1  <  IJ-tT^(t-u)F^i(X(u))dW(u)l^+1 
s  s 

+  f  I JUT  (u-v)F^^  (X(v)  )dW(v)  I  du 

J  J  a  -a+ 1  _ _ , 


1  /  y 

and  the  spatial  regularity  of  /  T  (•-s)F  (X(s) )dW (s) , as  proved  in  Kote- 

O 

lenez  [27]  (generalizing  a  result  of  Dawson  [ 10]  -  cf.  also  Kotelenez 
[25]).  That  the  spatial  regularity  in  (3.8)  is  maximal  follows  from  the 

f-  1/9 

Gaussianity  of  f  U  (t,s)F  .  (X(s))dW(s)  -  cf.  for  details  Kotelenez 
u  J  -a  -<x+ 1 
o 

[25]. 

(iii)  The  Markov  property  follows  from  Arnold,  Curtain  and  Kotelenez 
[  2  ],  (3.9)  follows  from  Curtain  [  6 ]. 

Remark  3 . 1 

1.  The  final  result  can  be  expressed  by 


(3.10)  XM  =  X  +  -  Y  +  o  ( - )  , 

‘  /vN 


where  X,  is  the  local  stochastic,  i.e.,  mezoscopic  description,  assuming 

X^  (O)  being  near  to  homogeneity  ,  X  is  the  deterministic  homogeneous 

state  solution  of  (2.1),  and  Y  is  the  mild  solution  of  (3.7)  which  is 

a  generalized  Gauss-Markov  process  if  Y  is  Gaussian,  and  o( -  )  is 

°  v'vN' 

the  error  term. 


IX.  Let  us  now  assume  that  we  do  not  start  in  (2.1)  with  a  constant  but 
with  some  other  positive  bounded  and  possibly  smooth  function  XD(q) . 


Then  the  difference  X  (t)  -  X(t)  satisfies 

N 

(3.11)  XM(t)  -  X{t)  =  F  (t)  +  G  (t) , 

N  N  N 

rt 

where  F^it)  is  the  r.h.s.  in  (2.13)  and  GN(t)  =  J  UN(t,s) (DA^-DA)X(s)ds. 

o 

Of  course,  (3.11)  will  also  tend  to  zero  under  the  assumptions  of  Theo¬ 
rem  3.1.  However,  in  view  of  (2.18)  we  must  normalize  (3.11)  by  multi- 

1/2 

plying  both  sides  by  (vN)  (modulus  a  constant)  in  order  to  obtain  a 

Gaussian  correction  term.  On  the  other  hand,  the  convergence  of 

(vN) 1  c_,(t)  to  zero  with  v  =  iP  and  p  >  1  does  not  hold(in  general)  in 

function  norms  and  for  p  <  2  Arnold  and  Theodosopulu  [  3 ]  have  shown 

in  the  one-dimensional  case  that  the  variance  of  Z  (t)  (the  martingale 

N 

part  of  X  ( t )  tends  to  <*>  in  L  -norm.  This  problem  and  related  questions 
N  2 

will  be  inves  igated  in  a  forthcoming  paper. 
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